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Abstract. A gauge-invariant formulation of constrained variational calculus, based on the intro- 
duction of the bundle of affine scalars over the configuration manifold, is presented. In the resulting 
setup, the "Lagrangian" ^ is replaced by a section of a suitable principal fibre bundle over the veloc- 
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Introduction 

Almost ten years ago, a mathematical setting for a formulation of Classical Mechanics, automat- 
ically embodying its gauge invariance, has been introduced 121 12- Besides that, a more recent paper 
[1] proposes a geometric revisitation of the calculus of variations in the presence of non-holonomic 
constraints. 

The present work uses the arguments of |1] and the geometrical framework provided by [21 [3] to 
analyze the underlying gauge structure of constrained variational calculus. As a by-product, the 
resulting scheme allows to obtain the relevant ingredients that are commonly used in the variational 
context as canonical geometrical objects. 

The topic will be developed within the family of differentiable curves alone, thus avoiding all the 
issues coming from the possible presence of corners which are of poor significance from the gauge- 
theoretical point of view. The extension to the piecewise differentiable case, surely interesting on 
its own, can however be easily pursued. 

Consider an abstract system B, subject to a set of differentiable and possibly non-holonomic 
constraints, and let A denote its admissible velocity space. 

Then, define an action functional by integrating a suitable differentiable "cost function" on A 
— called the Lagrangian — along the admissible evolutions of the system. In its essential features, 
the variational problem we shall deal with is the one of characterizing, among all those evolutions 
of B which fulfil the restrictions imposed by the assigned constraints and connect a fixed pair of 
configurations, the extremals (if at all) of the given action functional. As we shall see, the whole 
topic has very close links with optimal control theory. 

For the sake of convenience, the first part of the paper shall be devoted as a reference tool 
consisting of a brief review of a few basic aspects of jet-bundle geometry and non-holonomic 
geometry as well as of those contents of [H [21 [3] that will be involved in the subsequent discussion. 
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This will include, among other things, a revisitation of the Lagrangian and Hamiltonian bundles. 
All results will be stated without any proofs nor comments. 

Afterwards, we shall develop an algorithm able to establish a canonical correspondence between 
the input data of the problem, namely the constraints and the Lagrangian function, and an al- 
ternative free variational problem over a distinguished overlying affine bundle ^(^) — >■ In the 
resulting setup, under suitable hypotheses, the gauge-independent problem in '^^(A) is proved to 
be equivalent to the actual constrained one. This clarifies the geometrical essence of Pontryagin's 
method, based on the introduction of the costates. 

As the final element, the circumstances under which the correspondence between the solutions of 
the variational problem in A and those of its associated problem in '^{A) is 1 — 1 are investigated. 



1. Geometric setup 

1.1. Preliminaries. Let Vn+i — > R denote an (n+l)-dimensional fibre bundle, henceforth called 
the event space and referred to local fibred coordinates t, g^, ...,(?" . 

Every section 7: IR, ^ Vn+i, locally described as = q'''{t), will be interpreted as an evolution 
of an abstract system B, parameterized in terms of the independent variable t. The first jet- 
space ii(Vn+i) — ^ Vn+i is an affine bundle over Vn+i, modelled on the vertical space V{Vn+i)- 
By the very definition of jet-bundle, every section 7: IR — t- Vn+i may be lifted to a section 
ji(7): H ji{Vn+i), simply by assigning to each f G R the tangent vector to 7. The section 
ji(7) will be called the jet-extension of 7 on ji(Vn+i). 

Both spaces ji(V„,+i) and V{Vn+i) may be viewed as submanifolds of the tangent space T(Vn+i) 
according to the identifications 

(Lla) Ji(V„+i) = {ze T{Vn+i) I {z, dt) = 1} 

(1.1b) V{Vn+i) = {v G T(V„+i) I (v, dt) = 0} 

The terminology is borrowed from Classical Mechanics, where B is identified with a material 
system, the manifold Vn+i with its configuration space-time, the projection t : Vn+i — > R with 
the absolute time function and the jet-space ji(Vn+i) with the velocity space of B. 

The geometry of the manifold jiiVn+i) will be regarded as known. The reader is referred to 
[U H] for the notation, the terminology and a thorough analysis. Unless otherwise stated, given any 
set of local coordinates on Vn+i , the corresponding local jet-coordinate system on ji(Vn+i) will be 
denoted by t, g-^, . . . , g", g-^, . . . , g". 

The dual of the vertical bundle, henceforth denoted by V^*(Vn+i) — > Vn+i, is, in view of eg. (jl.lbp . 
canonically isomorphic to the quotient of the cotangent space T*(Vn+i) by the equivalence relation 



it{g) = vr(cr') 

G — g' (X dt|7r(o-) 



Every local coordinate system t,q^ in Vn+i induces fibred coordinates t,q^,pi in V*{Vn+i), with 
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and transformation laws 



(1.2) 



t = t + c, q' = q'it,q\...,q'') , Pi = Pk 



dq^ 
dq^ 



The annihilator of the tangent distribution to the totahty of the jet-extensions of sections 7 is 
a subspace C(ji(Vn+i)) of T*(ji(Vn+i)) , cahed the contact bundle. 

Alternatively, this last may be seen as the pull-back of the space V*{Vn+i) through the map 
ji(Vn+i) Vn+i- As such, C(ji(V„+i)) is, at the same time, a vector bundle over ji(Vn+i) and 
an affine bundle over V*{Vn+i) ■ The manifold C(ji(Vn+i)) will be referred to coordinates t, q', q^,Pi, 
related in an obvious way to those in ji(Vn+i) and in V*{Vn+i) ■ Every cr G C(ji(Vn+i)) will be 
called a contact 1-form over ji(Vn+i)- 

1.2. Non— holonomic constraints. Let A denote an embedded submanifold of ji{Vn+i), fibred 
over Vn+i ■ The situation, summarized into the following commutative diagram 



A 



Jl(Vn+l) 



(1.3) 



+1 



Vn 



+1 



provides the natural setting for the study of non-holonomic constraints. 

The manifold A is referred to local fibred coordinates t,q^ , . . . ,q"', , . . . , z'^ with transformation 
laws 



(1.4) 



t = t + c, r = q\t,q\...,q"), z^ = z^{t,q\ . . . ,q\ z\ . . . , z') 



while the imbedding i : A^ ii(Vn+i) is locally expressed as 
(1.5a) q^ = tl)\t,q^,...,q^,z^,...,z'^) , i = l,...,n , rank 

or, alternatively, may be implicit represented as 
(1.5b) g''{t,q\...,q^,q\...,q'')=0 , a = l,...,n-r, rank 



d{^'^ ■ ■ ■ ■!/>") 



9(^1 



= r 



d(g^ 



... g'. 



d{q 



n — r 



In the following, we shall not distinguish between the manifold A and its image i{A) C ji(V„+i). 

In the presence of non-holonomic constraints, an evolution 7: R — ?> Vn+i is called admissible 
if and only if its first jet-extension is contained in A, namely if there exists a section 7: IR ^ ^ 
satisfying ji(7r • 7) = i • 7. Expressing any section 7 in coordinates as = g*(t) , = z^{t) , the 
admissibility requirement takes the explicit form 



(1.6) 



^ = nt,q\t),...,q^{t),z\t),...,z'{t)) 



The concepts of vertical vector and contact 1-form are easily extended to the submanifold A: as 
usual, the vertical bundle V(A) is the kernel of the push-forward tt* : T{A) — >■ T(Vn+i) while 
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the contact bundle C{A) is the pull-back on A of the bundle C(ji(Vn+i)), as expressed by the 
commutative diagram 

C{A) > C(ii(V„+i)) > V*{Vn+i) 

(1.7) 

A — ^ ji(V„+i) — ^ V„,+i 
The latter allows to regard the contact bundle C{A) as a fibre bundle over the space V*{Vn+i), 
identical to the pull-back of V*(yn+i) through the map A Vn+i- 

1.3. Infinitesimal deformations of sections, (i) Quite generally, given a section 7 : IR, ^ Vn+i, 
a (weak) deformation of 7 is a 1-parameter family of sections 75 , C G {—e,e) depending differen- 
tiably on ^ and satisfying 70 = 7- 

In the presence of non-holonomic constraints, a deformation 7^ is called admissible if and only 
if each section 7^ : K, ^ Vn+i is admissible in the sense of ^1.21 In a similar way, a deformation 
7^ of an admissible section 7 : El ^ ^ is called admissible if and only if all sections 7^ : IR — ?> ^ 
are admissible. 

By definition, the admissible sections 7 : K, — Vn+i are in 1-1 correspondence with the admis- 
sible sections 7 : HI ^ ^ through the relations 

(1.8) 7 = 7r-7, ji{j)=i.j 

Every admissible deformation of 7 may therefore be expressed as 7^ = tt • 7^ , being 7^ : IR, — ?> ^ 
an admissible deformation of 7. 

In coordinates, preserving the representation 7 : = q^{t) , = z^{t) , the admissible defor- 
mations of 7 are described by equations of the form 

(1.9) 7^: q' = y^'{C,t), = CHi.t) 
subject to the conditions 

(1.10a) ^\Q,t)=q\t), C^{Q,t)=z\t) 

(1.10b) ^ = i;\t,ip\i,t)X\i,t)) 

For each t G IR, the curve ^ — ?• 7g(t) is called the orhit of the deformation 7^ through the point 
7(t) . The vector field along 7 tangent to the orbits at ^ = is called the infinitesimal deformation 
associated with 7^. 

Setting 



the infinitesimal deformation tangent to 7^ is described by the vector field 

while equation (jl.lObp is reflected into the relation 

dX^ d dip 



(1.13) 



dt dt 
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commonly referred to as the variational equation. 

The infinitesimal deformation tangent to the projection 7^ = tt • 7^ is similarly described by the 
field 

(1.") A- = ..*=f|^) m=X-wfA 



(ii) We now want to sketch out the construction of a suitable geometrical environment for the 
description of such infinitesimal deformations. In order to do so, we denote by V{^) B, the 
vector bundle over R formed by the totality of vertical vectors along 7, and by A (7) IR the 
analogous bundle formed by the totality of vectors along 7 annihilating the 1-form dt. Both bundles 
are referred to fibred coordinates — the former to t, and the latter to t, v^, — according to 
the prescriptions 



7(t(X)) 



XgA(7) ^ x = v\X)[^\ +w^{X)( ^ 



As proved in the first jet-bundle jiiVi'j)) is canonically isomorphic to the space of vectors 
along the jet-extension ^1(7) annihilating the 1-form dt. 

In jet-coordinates, the identification is expressed by the relation 

The push-forward of the imbedding A ^ ji(Vn+i), restricted to the subspace ^(7) C T{A), makes 
the latter into a subbundle of ji{V{^)). This gives rise to a fibred morphism 



(1.15a) 



^(7) Ji(^(7)) 



y(7) = y(7) 



expressed in coordinates as 



All previous results are then summarized into the following 

Proposition 1.1. Let 7 : R — t- Vn+i and 7 : IR ^ ^ denote two admissible sections, related by 
equation U.8\) . Then: 

i) the infinitesimal deformations of 7 and those of 7 are respectively expressed as sections 
X : E ^ V{j) and X -.R^ A{j); 

ii) a section X : R ^(7) represents an admissible infinitesimal deformation of 7 if and 
only if its first jet-extension factors through ^(7), i.e. if and only if there exists a section 
X : R ^ ^(7) satisfying ji{X) = i^X; conversely, a section X : R ^ ^(7) represents 
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an admissible infinitesimal deformation of 7 if and only if it projects into an admissible 
infinitesimal deformation of 7, i.e. if and only if i^X = ji(7r*X) . 

From a structural viewpoint, Proposition 11.1 1 points out the perfectly symmetric roles respectively 
played by diagram (jl.Sp in the study of the admissible evolutions and by diagram (jl.lSap in the 
study of the infinitesimal deformations, thus enforcing the idea that the second context is essentially 
a "linearized counterpart" of the former one. 

1.4. The gauge setup. 

1.4.1. The Lagrangian bundles. Given any system subject to (smooth) positional constraints, we 

7r t 

introduce a double fibration P — > V„+i — > R, where: 

i) Vn+i — > is the configuration space-time of the system; 
a) P —> Vn+i is a principal fibre bundle with structural group (R, +) . 
As a consequence of the stated definition, each fibre Px ■= 7r~^(x), x € Vn+i is an affine 
1-space. The total space P is therefore a trivial bundle, diffeomorphic in a non-canonical way 
to the Cartesian product Vn+i x IR, called the bundle of affine scalars over Vn+i- 

The action of (R, +) on P results into a 1-parameter group of diffeomorphisms P ^ P , 
conventionally expressed through the additive notation 

(1.16) ^p^{jy) := + ^ V^GR, z/GP 
Every map u : P ^ R satisfying the requirement 

u{u + = '^i^) + ? 

is called a (global) trivialization of P. If u, u' is any pair of trivializations, the difference u — u' 
is then (the pull-back of) a function over Vn+i- Moreover, every section q: Vn+i — > P determines 
a trivialization G ^[P) and conversely, being the relation between <j and expressed by the 
condition 

(1.17) V = ?(7r(i/)) + u^{v) Vi/ G P 

Therefore, once a (global) trivialization u: P ^ R has been chosen, every section <;^: Vn+i — )■ P is 
completely characterized by the knowledge of the function / = ?*(it) G ^iVn+i)- 

The assignment of u allows to lift every local coordinate system t,q^ , . . . over Vn+i to a 
corresponding fibred one t,q^, . . . ,q^,u over P, being the most general transformation between 
fibred coordinates of the form 

t = t + c, r = q\t,q\---,qn, u = u + f{t,q\ . . . ,q^) 

The action of the group (R, +) on the manifold P is expressed in fibred coordinates by the relations 

t{iy + i)=t{v), q\iy + i)=q'{v) , u{u + i) = u{v) + ^ 

As a result, the generator of the group action ()1.16p . usually referred to as the fundamental vector 
field of P, is canonically identified with the field ^ . 

The (pull-back of the) absolute time function determines a fibration P — ^ R whose associated 
first jet-space is indicated by ji(P, R) —> P and is referred to local jet-coordinates t,q\u,q^,u 
subject to transformation laws 

(1.18a) t = t + c, q' = q\t,q\...,q^) , u = u+ f{t,q\ . . . ,q^) 
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(1.18b) 



dq 



dq^ ~^ dt ' 



dq' ^ df df 

The manifold ji{P, R) is naturally embedded into the tangent space T{P) through the identifi- 
cation 

ji{P,R) = {zeT{P) I {z,dt) = 1} 
expressed in local coordinate as 



(1.19) 



z€ji{P,R) 



d ■ d ■ d 
- + ^X.)^ + n^(.)- 



■k{z) 



In addition to the jet attributes, the space ji(P, HI) inherits from P two distinguished actions of 
the group (R, +), related in a straightforward way to the identification (jl.l9p . 

The former is simply the push-forward of the action ()1.16p . restricted to the submanifold 
ji(P, R) C T{P). In jet-coordinates, a comparison with equation p.l9p provides the local repre- 
sentation 



d_ 

dt 



(l-20a) (^5)Jz) = - + 

expressed symbolically as 

(1.20b) (V^)^ : {t,q\uA\ul {t,q\u + ff^u') 

The quotient of ji(P, R) by this action is a (2n + 2)-dimensional manifold, denoted by £(Vn+i). 
The quotient map makes ji{P, R) into a principal fibre bundle over C{Vn+i) , with structural group 
(R, +). Furthermore, by equation ()1.20bp . C{Vn+i) is an affine fibre bundle over Vn+i with local 
coordinates t,q'^,q^,u. 

The latter action of (R, +) on ji{P, R) follows from the invariant character of the field ^ and 
is expressed in local coordinates by the addition 



(1.21a) 



d_ 

du 



d_ 

dt 



d 



d_ 

du 



7r(z) 



summarized into the symbolic relation 

(1.21b) cP^ : {t,q\u,q\e) {t,q\u,q\e + 

The quotient of ji{P, R) by this action is once again a (2n + 2)-dimensional manifold, denoted by 
C^iVn+i)- As before, equation ()1.21bp points out the nature of C^{Vn+i) as a fibre bundle over P 
(as well as on Vn+i), with coordinates t,q^,u,q^. The quotient map makes ji(P, R) — )■ C^{Vn+i) 
into a principal fibre bundle, with structural group (R, +) and group action (|1.21ap . 

Eventually, the group actions ()1.20ap . ()1.21ap do commute. Therefore, each of them may be 
used to induce a group action on the quotient space generated by the other one. This makes 
both C{Vn+i) and C^{Vn+i) into principal fibre bundles over a common "double quotient" space, 
canonically diffeomorphic to the velocity space ii(Vn+i). 

The situation is summarized into the commutative diagram 



+1) 



il(Vn+l) 
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in which all arrows denote principal fibrations, with structural groups isomorphic to (R, +) and 
group actions obtained in a straightforward way from equations (|1.20bp . (|1.21bp . The principal 
fibre bundles C{Vn+i) — > ji(V„+i) and C^iVn+i) — > Ji(V„+i) are respectively called the Lagrangian 
and the co~Lagrangian bundle over ji(Vn+i)- 

The advantage of this framework is most appreciated by giving up the traditional approach, 
based on the interpretation of the Lagrangian function as the representation of a 

(gauge-dependent) scalar field over jiiVn+i) and introducing instead the concept of Lagrangian 
section, meant as a section 



and so it does still rely on the assignment of a function ^{t,q^ ^q^) over ji(Vn+i)- However, as 
soon as the trivialization is changed into u = u + f , the representation (|1.22p undergoes the 
transformation law 



involving a different, gauge-equivalent, Lagrangian. 

1.4.2. The non-holonomic Lagrangian bundles. In the presence of non-holonomic constraints, the 
construction of the Lagrangian bundles may be easily adapted to the submanifold A, through a 
straightforward pull-back process. 

The situation is conveniently illustrated by means of a commutative diagram 




(L22) 



(L23) 



^ = ^ + / = if (t, q\ f) + / := if"(t, q\f) 





(L24) 



ii(v„+i) 





C{A) 



A 



where: 



C{A) and C'^{A) are respectively the pull-back of C{Vn+i) and £'^(Vn+i) on the submanifold 
A^ji{Vn+i); 
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• the space jj^(P, R) may be alternatively seen as the pull-back of the jet-bundle 
ji(P, R) — £(Vn+i) on the submanifold C{A) — )■ £(Vn+i) or as the pull-back of 
ji(P, R) ^ C%Vn+i) on C%A) ^ C%Vn+l). 

The geometrical properties of the above-defined pull-back bundles are straightforwardly inherited 
from their respective holonomic counterparts. In particular: 

• Every choice of a trivialization u of P allows to lift any coordinate system of A to coordi- 
nates t, q^, z^, u on C^iA) , t, g*, z^, u on C{A) and t, q\u, z^,u on j'f{P, R) . The resulting 
coordinate transformations are obtained by completing equations (|1.4p with (the significant 
part of) the system 



(1.25) u = u + f{t,q\...,q^) , i, = u + + ^^\t, q\ z"") := u + f 



• The embeddings C{A) CiVn+i) , C%A) /:^(V„+i) as weh as jf{P, R) ji{P, R) are 
all locally described by equation ()1.5ap . 

• Both actions (|1.2Uap . (|1.21ap of the group (R, +) on ji(P, R) preserve the submanifold 
ji'{P, R) thereby inducing two corresponding actions (^g)^ and 0g on j'fiP, R), expressed 
in coordinate as 

(1.26a) (V'^)^ : {t,q\u,z'^,u) {t,q\u + ^, z^,u) 

(1.26b) (p^ : {t,q\u,z^,u) {t,q\u, z'^,u + 

Acting in the same way as before, it is easily seen that the manifold ji'{P, R) is a principal 
fibre bundle over C{A) under the action (V'c)^, ■, as well as a principal fibre bundle over C'^{A) 
under the action (f)^. Moreover, both C{A) and C^iA) are principal fibre bundles over A 
under the (induced) actions (V'^)^ and (p^ respectively. Accordingly, all arrows in the front 
and rear faces of diagram ()1.24p express principal fibrations, while those in the left and 
right-hand faces are principal bundle homomorphisms. 

Preserving the terminology, the principal bundles C{A) — )• A and C^iA) — >■ A are respectively called 
the non-holonomic Lagrangian bundle and the non-holonomic co-Lagrangian bundle over A, while 
any section £: A ^{A) is referred to as a non-holonomic Lagrangian section. 

Once a trivialization u of P has been fixed, the description of £ takes the local form 

(1.27) u = .^{t,q\z^) 
which undergoes the transformation law 

(1.28) i, = u + j = ^(t, q\z^) + ^ + ^^^i := ^'(t, q\ z^) 
under an arbitrary change u — ?> u + f{t,q^ , . . . ,q^). 

1.4.3. The Hamiltonian bundles. Parallelling the discussion in S I1.4.H we shall now briefly go over 
the construction of the Hamiltonian bundles on Vn+i • To this end, we focus on the fibration 
P — > Vn+i, and denote by vr: ji(P, Vn+i) — > P the associated first jet-space. 
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Every fibred coordinate system t,q^,u on P induces local coordinates t,q^,u,po,pi on ji{P, Vn+i), 
with transformation group 

(1.29a) t = t + c, r = q'{t,q\...,q^), u = u + f{t,q\ . . . ,q^) 

df / df\ dq^ ( df\ dq^ 

(1.29b) p-o=Po + ^+^P. + ^J ^, m=[p, + ^,) g- 

The manifold ji{P,Vn+i) is naturally imbedded into the cotangent space T*{P) through the 
identification 

ii(p,v„+i) = lr?Gr*(p) / ^ 



du / 

expressed in local coordinate as 

(1.30) 77 e ji(P, V„+i) ^ T] = [du- Po{v)dt - pMdq']^^^^ 

Furthermore, the jet-bundle structure endows ji{P, Vn+i) with a contact bundle, locally generated 
by the Liouville 1-form 

(1.31) e = du - podt - pidq' 

which is an intrinsically defined object, as ensured by equations ()1.29l a.b). 

Exactly as in the Lagrangian case, one can easily establish two distinguished actions of the group 
(IR, +) on ji(P, Vn+i), expressed locally as 

(1.32a) ii'dAv) ■■= ii^-O^iv) = [du-po{v)dt-p,{rj)dq']^^^^^^ 

(1.32b) (/>5(r/) := rj-C{dt)^i^) = [du - {po{7]) + dt - pi{7])dq']^^^^ 

In this connection, we point out that: 

• The direct product of the actions ()1.32al b) makes ji{P, Vn+i) into a principal fibre bundle 
over a (2n + l)-dimensional base space n(V„+i), with coordinates t,q'^,pi, called the phase 
space. 

• In view of equations ()1.2p . ()1.29l a.b). the phase space n(Vn+i) is readily seen as an affine 
bundle over Vn+i, modelled on V*{Vn+i)- 

• The quotient of ji(P, Vn+i) by the action ()1.32ap . denoted by ?^(Vn+i), is an affine bundle 
over Vn+i, modelled on the cotangent space r*(V„+i) and called the Hamiltonian bundle. 

• Any trivialization n : P — )• R allows to lift every local coordinate system t,q^, . . . ,q" on Vn+i 
to a corresponding one t, q^ , . . . ,q"^,po,pi, . . . ,pn on ?^(Vn+i) , subject to the transformation 
law 

(1.33) = + _ , p.=p. + _. 

further to a change of u into u = u + f{t, q^, . . . , q^) . 

• The quotient map makes ji(P, Vn+i) into a principal fibre bundle over ?^(Vn+i), with 
structural group (R, +) and fundamental vector ^ . 

• The canonical Liouville 1-form ()1.3ip endows ji(P, Vn+i) ^(Vn+i) with a distinguished 
connection, called the Liouville connection. 
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• The action ()1.32bp "passes to the quotient", thereby making TiiVn+i) into a principal fibre 
bundle over the phase space n(Vn+i)- 

• The quotient of ji(P, Vn+i) by the action ()1.32bp . denoted by 'W^iVn+i), is a 
(2n + 2)-dimensional manifold, with coordinates t,q'^,u,pi, called thek co-Hamiltonian 
bundle. The quotient map makes ji{P,Vn+i) into a principal fibre bundle over 'H'^{Vn+i)- 

• The action ()1.32ap . suitably transferred to T-L^iVn+i), makes the latter into a principal fibre 
bundle over n(Vn+i) . 

The previous discussion is summarized into the commutative diagram 

ji(p,v„+i) > n^Vn+i) 

(1.34) 

nvn+i) > n(v„+i) 

all arrows denoting principal fibrations with structural group isomorphic to E,. 

As implicit in the notation, the manifold ji {P, Vn+i) is indeed identical to the pull-back of T-L'^{Vn+i) 

over T-l{Vn+i), as well as the pull-back of TiiVn+i) over 'H'^{Vn+i)- 

1.5. Further developments, (i) The identifications (jl.l9p . (jl.30p provide a natural pairing be- 
tween the fibres of the first jet-spaces ji(P, R) P and ji(P, Vn+i) — ^ P, locally expressed 
as 



^1.35) {z,v) 



d d d 

- + <^^(.)^ + n(.)- 



, [du - po{r])dt - pi{'n)dq']^^^^ 



for all z G ji{P, R) , r/ G ji{P, V„+i) satisfying Tr{z) = 7r(r/) . 

In view of equations (|1.20ap , (|1.32ap , the correspondence (|1.35p satisfies the invariance property 

(1-36) (ii^Mz) , ii^Mv)) = 

thereby inducing an analogous pairing operation between the fibres of the bundles £(Vn+i) — > Vn+i 
and 7i{Vn+i) Vn+i, or — just the same — giving rise to a bi-affine map of the fibred product 
£(Vn+i) xv„+i ^(Vn+i) onto R, expressed in coordinates as 

(1.37) Cfi FiCfi) := u{0-po{fi)-p,{fi)q\0 

(ii) Let S denote the submanifold of C{Vn+i) xv„+i T~{-{Vn+i) described by the equation 
(1-38) 5 = {(C,^) e £{Vn+i) xv„+i n{Vn+i) I FiCfi) = O} 

A straightforward argument, based on equation ()1.37p . shows that the submanifold S is at the 
same time a fibre bundle over CiVn+i) as well as over ?^(Vn+i). The former case is made explicit 
by referring S to local coordinates t,q^ ,u,pi, the pis been regarded as fibre coordinates. The 
latter circumstance is instead accounted for by referring S to coordinates t,q^,q^,po,pi, related to 
the previous ones by the transformation 

ii = po + Piq' 
and with the g*'s now playing the role of fibre coordinates. 
The restriction to the submanifold S of the action 

<p^iC,fi) := {<p^iO,Mi^)) v(c,/i) G ciVn+i) xv„+, niVn+i) 
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makes the latter into a principal fibre bundle over a (3n + 1) -dimensional base space '^(ji(Vn+i)), 
with coordinates t,q'^ ,pi. 

Depending on the choice made for the local coordinates over S, the resulting group action may 
be expressed symbolically either as 



(1.39a) 
or 

(1.39b) 



: {t,q\q\u,pi) — > {t,q\q\u + i,pi) 



The situation is summarized into the following diagram 



(1.40) 



-HiVn 



+1; 



n(v. 



n+lj 



Ji(V, 



n+1) 



Vn+1 



In view of equations (|1.2|) and (|1.29l b). the manifold ^{ji{Vn+i)) is — by construction — an affine 
bundle over ji(Vn+i), modelled on the contact space C{ji{Vn+i)) ■ 



(ill) The canonical contact 1-form ()1.3ip of ji{P, Vn+i) can be pulled-back onto the fibred prod- 
uct ji(P, IR) Xp ji(P, Vn+i) , thus endowing the principal fibre bundle ji(P, R) Xp ji{P,Vn+i) 
C{Vn+i) Xv„+i ^(Vn+i) with a canonical connection. 

For every choice of the trivialization u of P Vn+i , the difference du — Q is (the pull-back of) 
a 1-form G)„ on CiVn+i) Xv„+i 'HiVn+i) , locally expressed as 

(1.41) G„ = podt + pidq' 
and subject to the transformation law 

(1.42) = (^po + dt+ (^p, + dq' = Qu + df 

under an arbitrary change u — > n = n + /(t, q^, . . . , g") . 

The form 0^ can now be once again pulled-back onto S . In this last step, depending on the 
choice of the coordinates over S , the resulting 1-form is locally expressed as 



(1.43) 



©M = Podt + Pi dq^ = iidt + Pi [dq^ — q^ dt) 



Hence, the submanifold S is provided with a distinguished 1-form , defined up to the choice of 
the trivialization of P . 
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(iv) In the presence of non-holonomic constraints, the left-hand face of ()1.40p , 
(1.44) 



'^(jl(Vn+l)) > Jl(Vn+l) 

may be easily pulled-back through the imbedding A Ji(Vn+i), giving rise to the analogous 
diagram 

(1.45) 

^(^) > A 

By construction, the manifold 5"^ is then a principal fibre bundle over the base space 'W{A) under 
the (induced) action 

(1.46) (j)^ : {t,q\z'^,u,pi) {t,q\ z'^,u + ^Pi) 

while, in the same manner as before, 'if (A) is an affine bundle over A modelled on the non- 
holonomic contact bundle C{A). 

By means of the pull-back procedure, the canonical form (|1.43p determines a distinguished 
1-form on 5"^, locally expressed by0 

(1.47) 9„ = podt+pidq' = iidt+pi {dq' - ip' dt) 

(v) Every non-holonomic Lagrangian section i: A determines a trivialization 

ip£: C{A) ^ R of the bundle C{A) A. Let (p£ := ^^{(pi) denote the pull-back of ^pi to 5"^, 
locally expressed as 

(1.48) ^i{t,q\z^,u,pi) = ipe{t,q\z'^,u) = ii- ^{t,q\z^) 

From this, taking equation (|1.46p into account, it is an easy matter to check that the function pi is 
a trivialization of the bundle <S"^ — )• ^{A) and that, as such, it determines a section i: 'i^(A) 5"^, 
locally described by the equation 

(1.49) u = ^{t,q\z^) 

In brief, every section £: A ^ ^(^) be lifted to a section £: 'ia{A) — ?> S-^ . The local 
representations of both sections are formally identical and they obey the transformation law ()1.28p 
for an arbitrary change of the trivialization u: P ^ K, . 

The section £ : ^{A) — > S-^ may now be used to pull-back the form (jl.47p onto ^{A) , hereby 
getting the 1-form 

(1.50) Gppc := r(G„) = ^dt + pi {dq' - ip' dt) := -J^ dt + pidq' 

henceforth referred to as the Pontryagin-Poincare-Cartan form. 

Needless to say, the difference := Piip^ — , known in the literature as the Pontryagin 
Hamiltonian, is not an Hamiltonian in the traditional sense. 



-'^The same symbol 0„ will stand for both the form (|1.43p and its pull-back on A. 
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Remark 1.1. The nature of the Pontryagin Hamiltonian may be understood by pointing up that, 
in view of equations ()1.5p . ()1.37p . ()1.38p . the space 5"^ is — by construction — a submanifold of 
jC(Vn+i) xv„+i T~l-{Vn+i) locally described by the equations 

F(C,/i) = 5"(t,'?(C),'7(C)) = , a = l,...,n-r 

for any {(, fi) G C{Vn+i) Xv„+i TiiVn+i) ■ 

Hence, the manifold S-^ may be equivalently referred to both local coordinates t,q^,pi,z^,u and 
t,q^,Pi,z^,Po, related one another by the transformation 

ii = Po + Pi'tp'{t,q^,...,q'',z^,...,z'') 

In the former circumstance, the section £: 'io'{A) — > S'^ is locally represented by equation ()1.49p . 
while in the latter case its local representation involves the Pontryagin Hamiltonian Jf{t,q\pi,z^) 
in the form 

(1.51) PO = -Jeit,q\p^,z^) 

2. Application to Constrained Variational Calculus 

2.1. Problem statement. As already mentioned, we shall consider a constrained abstract system 
B. Given a differentiable function _Sf € ^{A) on the space A (called the Lagrangian) and denoted 
by 7 the lift to A of an admissible evolution 7 of the system, define the action functional 

(2.1) Z[7] := [^{t,q\...,q",z\...,z')dt 

We intend to characterize the extremals ( if at all ) of the functional I [7] among all the admissible 
evolutions 7: [toj^i] — ^ Vn+i connecting two given configurations. 

As we shall see, it turns out to be more convenient to approach the topic by defining an alternative 
variational problem on the manifold ^{A) and proving that, under suitable hypotheses, the latter 
is equivalent to the original one. 

Remark 2.1. The constrained variational problem based on the functional i2. 1\) may be also viewed 
as a typical optimal control problem. 

As in §[23) ths admissibility of a given evolution of B is expressed in local coordinates by equation 
il.6\) . expressing the derivatives ^ in terms of a smaller number of variables z"^, A = 1, . . . ,r. 

Therefore, every concurrent assignment of both the values of z^{t), . . . , z^{t) and of a point in the 
event space Vn+i determines an admissible evolution of the system as the solution of the ordinary 
differential equations U.6\) with the given initial conditions. This makes the z^ 's into the controllers 
of the evolution and, as such, they usually go in the literature under the name of controls. 

In this sense, the search for the curves 7 = 7(t) along which the functional h2. i|) takes its extremal 
values may be equivalently seen as the one for those particular controls which optimize the evolution 
of the system. 

Actually, in the absence of specific assumptions on the nature of the manifold A, the functions 
z'^{t), in themselves, have no invariant geometrical meaning. In this respect, attention should rather 
be shifted on sections a : Vn+i A , locally expressed as z"^ = z'^{t,q^, . . . , g") . 
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2.1.1. The gauge structure. Given any pair of 1-forms ^dt and J^' dt over A, their respective 
action integrals ^[7] = J^^dt and I' [7] = J^^^'dt give rise to the same extremal curves if the 
difference — ^) dt is an exact differential. 

This is easily seen as, under this circumstance, the equality ^ ^dt = § J^'dt holds along any 
closed curve, thereby entailing the relation 

^'[75] -Ali] = I {^'-^)dt^ [ {^'-^)dt 

•^75 -^7 

for any deformation 7^ vanishing at the end-points, whence also 

^(x'[7d -:i[7d) =0 

As a consequence, as far as a variational problem based on the functional (j2.ip is concerned, the 
Lagrangian function ^ G ,!^{A) is defined up to an equivalence relation of the form 

(2.2) ^ ^'_^ = ^, /g ^(Vn+l) 

Otherwise stated, the real information is not brought so much by ^ in itself as by a whole family 
of Lagrangians, equivalent to each other in the sense expressed by equation ()2.2p . 

The significance of the arguments developed so far relies actually on the fact, explicitly pointed 
out by equations (|1.27|) . (|1.28p . that the representation of an arbitrary section i: A ^ ^{^) 
involves exactly this family of Lagrangians, henceforth denoted by A{£) . A straightforward check 
shows that a necessary and sufficient condition for two sections i and i' to fulfil A(^) = A(£') is 
that the difference £' — i, viewed as a function over A , be itself of the form 

(2.3) r -I = fe ^{Vn+i) 

Thus we see that, within our geometrical framework, the equivalence relation (|2.2p between func- 
tions is replaced by the almost identical relation (j2.3p between sections. Intuitively, the latter is 
a sort of "active counterpart" of the transformation law (|1.28p for the representation of a given 
section l under arbitrary changes of the trivialization u: P ^ H. 

2.2. Extremals. To start with, we observe that the algorithm described in SI 1.5 1 allows to deduce 
— in a canonical way — a Pontryagin-Poincare-Cartan form 

(2.4) eppc = .^dt+p, {dq' - ip' dt) = -Jifdt + p., dq' 

on the manifold ^{A) from the only knowledge of the input data of the problem, namely 

i) the non-holonomic constraints, described by the imbedding i: A ^ Ji(Vn+i) and locally 
expressed by the equations (|1.5ap ; 

ii) the non-holonomic Lagrangian section £, locally represented in the form ()1.27p . being 
J^{t,q^ , z"^) the given Lagrangian function. 

To understand the role of the form (|2.4p in the present constrained variational context, we next 
focus on the fibration ^{A) — ^ Vn+i, given by the composite map '^{A) n(V„+i) — > Vn+i 
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and we define an action integral over ^(A), assigning to eacli section 7 : Vn+i — ^ '^(A), locally 
expressed by = q^t), = z^{t),pi = Pi{t) , the real number 



(2.5) Ml] ■■= j Qppc = [pi^-M'jdt 

Given any deformation 7^ preserving the end-points of f ■ 7 and indicating with 

X = X*(t)(^)- + T'^it) + ni(t)(^)- the corresponding infinitesimal deformation, we 
get the relation 

dX[%] 



5=0 J to 



L [\dt dp J''' 5g^y 



dt 



From the latter, taking the conditions X^{to) = X'^{ti) = into accomit, it is easy to conclude 
that the vanishing of ^ | under arbitrary deformations of the given class is mathematically 
equivalent to the system 

(2.6.) 1= a^=^ (^'^'^ ) 

^ dje _ dtp' _ Q 

dz"^ dz"^ dz^ 

Equation (|2.6a|) shows that the extremal curves of the functional (|2.5p are kinematically admissible. 
As they are extremals with respect to arbitrary deformations vanishing at the end-points, this 
automatically makes them extremals with respect to the narrower class of admissible deformations 
as well. 

Therefore, we can state that every "free" extremal of the functional ()2.5p projects onto an extremal 
curve 7:5* = q^{t) of the original problem. 

Conversely, given any admissible solution 7 of the assigned variational problem, it seems be- 
forehand hard to establish if and under which hypotheses there exists at least one extremal 7 of 
the functional (j2.5|) projecting onto 7. Heuristically, the associated variational problem in ^{A) 
may be viewed as the study of the functional ()2.ip where the kinematical admissibility condition 
p.6p does not play the role of an a priori request upon sections anymore but is instead retrieved 
afterwards by the method of Lagrange multipliers. As a consequence, it can be reasonable that, 
under suitable hypotheses, the equivalence between the two variational problems could be proved. 
Let us investigate this point. 

Referring to [1] for details, we recall that an admissible section 7: [toj^i] ^ Vn+i is called ordi- 
nary if and only if all its infinitesimal deformations vanishing at the end-points are tangent to a 
finite deformation with fixed end-points. A remarkable subclass of the ordinary sections is formed 
by the normal ones. We will dedicate the final paragraph to discuss their role in the present context. 



We now state the following 

Theorem 2.1. Every ordinary extremal 7 of the functional i2.1\) is the projection of at least one 
extremal 7 of the functional \2. 5|) . 
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Proof. Given an ordinary extremal 7 of (|2.ip . let 7 : = Q^{t), = z^{t) be its lift into 

A. As discussed in ^1.31 . an admissible infinitesimal deformation of 7 is then a vector field 
X = X^{t) {-^)xf + r'^(t) satisfying the variational equation 



dt \dq^)x^ \dz^J- 



more suitably written as 

(2.8) i.(A;x') = A',l^r- 

being A'^j{t) any non-singular solution of the matrix equation 
Furthermore, setting X*(to) = 0, the validity of the relation 

(2.0, (Ili)^r-^.. 

is then a necessary and sufficient condition for X^(ti) to vanish as well. The whole class of the 
infinitesimal deformations of 7 vanishing at its end-points is therefore in bijective correspondence 
with the totality of vertical vector fields r'^(t) {^x) - defined along 7 and satisfying the condition 

We next introduce n new functions pi = Pi{t) subject to the conditions 
(2-11) -TT + Pfc 



dt \ dq^ / - \ dq^ 

In view of equation ()2.9p . the functions Pi{t) are defined up to a transformation 
(2.12) pi{t) pi{t) = piit) + pjAh{t) 

with (/?!, . . . ,/3„) € R". Therefore, the thesis is proved as soon as we show that the stated hy- 
potheses entail the existence of at least one choice of (/3i, . . . , such that the resulting functions 
Piit) could be used to lift the curve 7 to an extremal 7: [to, ti] — )• 'ia[A) of the functional (|2.5p . 



(2.13) 



5=0 Jto \ dq^ dz^)^ 



Now, the extremality of the ordinary evolution 7 means that 

d^ 

for every infinitesimal deformation of 7 vanishing at its end-points. 

Moreover, taking the conditions X*(to) = X^{ti) = as well as the equations (|2.7p into account, 
we have 
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Substituting into equation ()2.13p , we conclude that 



(2.14) 



to 



Pk 



dip'' 



+ 



dt 







for each = T^{t) fulfilhng the condition (|2.10|) . In order for this to happen it is then necessary 
and sufficient the validity of the linear relation 



-Pk 



+ 



f^kA' 



dz^ 



7 \ dz^ J ^ \"~ / y 

or — what is the same — the existence of at least one solution Pi{t) = Pi{t) + PjA^i of the system 
(I2TT]1 such that 



(2.15) 



Pk 



dz^ 







Equations ()2.1ip and (|2.15p . together with the kinematical admissibility condition (true by hypoth- 
esis) are exactly the Euler-Lagrange equations for the functional (j2.5p . Therefore, the resulting 
curve 7:9* = 9*(i), z"^ = •z^(i); Pi = Pi{t) is an extremal of the functional ()2.5p which projects 
onto 7. □ 

Collecting all results, we have just shown that, as far as the ordinary evolutions are concerned, the 
original constrained variational problem in the event space is equivalent to a canonically associated 
free one in 'io{A) . 



Remark 2.2 (Same problem, equivalent solution). There is another possible approach to the prob- 
lem outlined in ^2.1\ which is slightly different but completely equivalent to the one taken so far. 
Apparently it just complicates matters without giving any significant advantage. On the other hand, 
it seems to be the most faithful translation of the original Pontryagin's treatment of the topic |6] 
into the present geometrical context. Hence — at least for historical reasons — it is worth telling 
about. 

The algorithm relies on the following considerations: 

(i) Given any section 7: [toi^i] ~^ P of the bundle of affine scalars, let 7 denote the restriction 
to C{Vn+i) of its jet-extension. 

The input data of the assigned problem are taken into account by the introduction of a notion of 
admissibility for 7 . This is accomplished by requiring the jet-extension 7 to belong to a submanifold 
A of C{Vn+i) , locally described by the equations 

(2.16) q' = i;\t,q\z^) , u = ^{t,q\z^) 

In other words, the simultaneous assignment of both the kinetic constraints and the Lagrangian 
function are used to express the submanifold A as the image i{A) . 

In this way, every admissible section q^ = q^{t) in Vn+i determines, up to the initial value u(to); 
a corresponding admissible section q^ = q^{t), u = u{t) of P . 

Compared with the main approach, the present formulation just replaces the section A ^ ^(-^) 
with the image space A = i{A) , viewed as a submanifold of C{Vn+i) ■ As we have seen, the 
submanifold A and, consequently, the section I are regarded as data of the problem and therefore 
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the representation u = ^{t,q^ , z^) is affected only by passive gauge transformations. The same 
variational problem is therefore bred by different submanifolds related one another by the action of 
the gauge group. 

(a) Let us consider the constrained variational problem with fixed end-points on the manifold 
C{Vn+i) , based on the functional 

(2.17) X[7] := I iidt 



7 denoting the jet-extension of an admissible section 7: [to, ii] — > P- The stated problem does not 
depend on a particular choice of the gauge, as the 1-form iidt is well-defined in £(Vn+i) up to a 
term fdt. 

In local coordinates, setting 7: = , u = u{t) , we have 

ii dt = u{ti) — u{to) 



and so, being the values of q^{to) and q^{ti) already fixed by the boundary conditions, the problem 
consists in finding a curve 7 over which the increment u(ti) — u{to) becomes stationary and whose 
projection onto V„+i joins the assigned end-points. 

(in) The submanifold A — )• C(yn+i) is lifted up onto a submanifold ^{A) — )• S whether by 
identifying '^(A) with the image space i{'if{A)) or by pulling back S onto A by means of the 
commutative diagram 

'r(i) — ^ S 

(2.18) 

A C{Vn+l) 

All the same, the embedding 'if (A) S is fibred onto n(Vn+i) and its expression in coordinate 
is formally identical to equations (|2.16p which are involved in the representation of the submanifold 
A. 

It is now possible to make use of the form ()1.43p to provide the manifold '^{A) with the 1-form 

(2.19) J* (9 J = ^dt + pi {dq' - dt) 

and, as a consequence, to define an action integral by the integration of (|2.19|) along any section 
7: [to,ti] ^^(i). 

Once again, this merely reproduces in the image space ^(^) the construction previously carried 
on in ^ \1.5[ In other words, the 1-form (|2.19|) is simply the image of the Pontryagin-Poincare- 
Cartan form ()1.50p under the diffeomorphism I: '^(A) — > ^{A). 

2.3. The case of normal sections. A significant role in the study of the variational problem 
based on the functional ()2.5p is played by the choice of the Lagrangian section as it = f{t,q), which 
is related to the particular situation determined by the ansatz ^{t,q\z^) = 0. 
Under the stated circumstance, the functional 

(2.20) Zo[7] := (^^'^ - ^"^dt 
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may therefore be used to single out a purely geometrical variational problem in the manifold 'if (A) . 
The corresponding extremal curves are easily seen to satisfy the Euler-Lagrange equations 

^-^•<'-'-^). . 

In view of equation (jl.6|) . every admissible evolution 7: [^Oi^i] ^ V^+i is the projection of at 
least one extremal 7 of the functional ()2.20p . For this reason, the projection of 7 under the map 
'j^{A) A coincides with the lift 7 : [toj^i] — ^ A. 

The extremals of ()2.20p are therefore in a bijective correspondence with the solutions Pi{t) of 
the homogeneous system ()2.2ip . with the functions q^t), regarded as given. In other words, 

the totality of extremals 7 of the functional (|2.2Up form a finite dimensional vector space over H, 
whose dimension will be referred to as the abnormality index of 7 := u • 7. 

Definition 2.1. An admissible curve 7: [toi^i] ~^ Vn+i is called normal whenever its abnormality 
index vanishes. Otherwise 7 is said abnormal. 

As a consequence, a normal curve 7: [^Oj^i] ~^ Vn+i is the projection of a unique extremal of 
the functional (|2.20p . namely of the curve 

7 : = q\t) , = z^it) , p,{t) = 



Coming back to the study of the main variational problem based on the action integral (|2.ip , we 
can now state 

Theorem 2.2. Every normal extremal 7 of the functional 12. 1\) is the projection of exactly one 
extremal 7 of the functional \2. 5|] . 



Proof. The crucial point, that we will not demonstrate here, is that the normal curves form a subset 
of the ordinary onesl^ 

In view of Theorem 12.11 this entails the extremal 7 to be the projection of at least one extremal 
7: = z^ = z^{t)^ Pi = Pi{t) of the functional (|2.5p . We still have to prove its uniqueness. 

To this end, we suppose the existence of a second extremal projecting onto 7, expressed in coordinate 
as 7': = q\t), z^ = z^{t), Pi = Ti{t). 

Taking equations (j2.2ip into account, it is easily seen that the contemporaneous validity of the 
Euler-Lagrange equations for both the curves 7 and 7' 

dq^ ,■ A, dpi di;^ d££ dip'' 

dq' A, dn d^'' 5if dip'' 

makes the curve = z^ = z^(t), pi = pi{t) — Ti{t) an extremal of the functional ()2.20p . 
Since 7 is — by hypothesis — a normal curve, this implies Pi{t) = Ti(t) and hence 7 = 7'. □ 



^See [T], Appendix B. 
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